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Abstract In these lecture notes we derive a generic holographic string theory real- 
ization of a p-wave superconductor and superfluid. For this purpose we also review 
basic D-brane physics, gauge/gravity methods at finite temperature, key concepts of 
superconductivity and recent progress in distinct realizations of holographic super- 
conductors and superfluids. Then we focus on a D3/D7-brane construction yielding 
a superconducting or superfluid vector-condensate. The corresponding gauge theory 
is 3+1 -dimensional jY = 2 supersymmetric Yang-Mills theory with SU(N C ) color 
and St/ (2) flavor symmetry. It shows a second order phase transition to a phase in 
which a U(l) subgroup of the SU(2) symmetry is spontaneously broken and typ- 
ical superconductivity signatures emerge, such as a conductivity (pseudo-)gap and 
the Meissner-Ochsenfeld effect. Condensates of this nature are comparable to those 
recently found experimentally in p-wave superconductors such as a ruthenate com- 
pound. A string picture of the pairing mechanism and condensation is given using 
the exact knowledge of the corresponding field theory degrees of freedom. 
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1 Introduction 

In these lecture note£] we genetically construct a holographic p-wave superconduc- 
tor. This introductory section serves to explain why this particular setup is of interest 
both from a string-theoretical as well as from a condensed matter physics point of 
view. In order to guide the unexperienced reader in section [2] we will describe the 
basic concepts in words and motivate the project. Detailing this overview we build 
the full holographic setup from scratch in a self-contained manner in section|3] Sec- 
tion |4] briefly introduces common holographic methods and then summarizes the 
major results known for the thermodynamics and the spectrum of D-brane systems. 



1 These lecture notes are based on joint research work with Martin Ammon, Johanna Erdmenger, 
Patrick Kerner and Felix Rust. 
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With these outcomes in the back of our minds we will understand how our system 
develops a new ground state and why it can be interpreted as a p meson superfluid 
or superconductor. We compute and discuss thermodynamic observables, conduc- 
tivities and the spectrum of our holographic superconductor in section [5] Finally, 
we put together all of our findings in order to draw a string theory picture for the 
pairing mechanism in our p-wave superconductor. These notes are self-contained 
and widely complementary to the notes focusing on holographic s-wave supercon- 
ductors CD 120. 



1.1 String motivation 

The setup of intersecting D-branes and especially the particular D3-D7 construc- 
tion presented here is interesting because of its diverse applications. It has been 
successfully used to model strongly coupled particle physics phenomena such as 
a deconfinement, or rather meson melting phase transition for fundamental matter 
and transport coefficients in the quark gluon plasma experimentally created at the 
RHIC collider in Brookhaven (see |4| for an introductory review and further ref- 
erences). On the other hand it has recently been used to model strongly correlated 
electron systems as they appear in superfiuids and superconductors in the realm of 
condensed matter physics |5]|6]. All these applications are also crucial checks of the 
basic principles and methods coming from the conjectured gauge/gravity correspon- 
dence. Receiving physically meaningful outcomes when applying these holographic 
methods to various systems accessible by experiment, strengthens our confidence in 
the gauge/gravity conjecture. 

From the string-theoretic point of view this particular setup is interesting be- 
cause of its simplicity, uniqueness (explained below) and the naturalness with which 
the symmetry is broken spontaneously here. The flavoid superfluid/superconductoiQ 
examined here was the first generic top-down string theory realization of a su- 
perfiuid/superconducting phase. In contrast to this the pioneering papers on holo- 
graphic superconductors Q E HI had exclusively treated gravity toy models which 
were not directly obtained from string theory, i. e. bottom-up models. Therefore the 
gauge/gravity correspondence could not be used to identify the exact gauge theory 
dual. This fact obstructed the interpretation of the outcomes. Furthermore these toy 
models had few restrictions on their parameters such that in principle a large pa- 
rameter space needed to be scanned, see e. g. iflOl . Our string-derived flavor super- 
fluid/superconductor overcomes those two problems: First the field theory degrees 



2 The term "flavor" superconductor stems from earlier applications of this setup to model strongly 
correlated high energy systems such as the quark gluon plasma. In the present case the name is 
not important and possibly misguiding, since it is really only essential that the system has a non- 
Abelian SU(2) symmetry. 

3 We are using the terms superfluid and superconductor interchangably here. For the considered 
phenomena this distinction does not make any difference. Some details to this distinction are given 
in section [2TI 
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of freedom are exactly known since it is simply jY = 2 supersymmetric Yang-Mills 
theory with an SU (2) flavor-symmetry. Second, the values of parameters, such as for 
example the dimension of the condensing operator, in our setup are severly restricted 
by their string theoretic derivation. In this sense this setup is "unique" compared to 
the big parameter space to be scanned in bottom-up toy models. Later, other top- 
down string realizations have been suggested and for example involve a consistent 
truncation of type IIB supergravity with a chemical potential for the R-charge ifTTl . 
and domain-wall solutions interpolating between AdS solutions with distinct radii 
which may be lifted to IIB supergravity or eleven-dimensional supergravity lfl2l . 



1.2 Condensed matter motivation 

From a condensed-matter physics point of view our flavor superconductor is highly 
interesting because it reproduces features which have been measured in experiments 
with unconventional superconductors, such as p-wave superconductivity, a system 
of strongly-correlated particles, a pseudo-gap in the frequency-dependent conduc- 
tivity (found in high temperature d-wave superconductors). Most of these phenom- 
ena lack a widely-accepted microscopic explanation by conventional approaches, so 
there is the hope that gauge/gravity can shed some light on the nature of these sys- 
tems. These systems usually contain strongly correlated electrons, so the dual weak 
gravity description is in principle accessible. Our system has a vector operator which 
condenses upon breaking a residual Abelian flavor symmetry spontaneously. This 
gives the vector order parameter as described below. So there is a preferred spatial 
direction in the superfluid/superconducting condensate. This is exactly the situa- 
tion recently found experimentally in the p-wave (explained below in section 12. Il l 
superconductor Sr2RuC>4 |13]. These materials are investigated with great excite- 
ment in the condensed-matter community because they are hoped to be usable for 
quantum computing fT4l . The reason is that the p-wave structure implies the pres- 
ence of non- Abelian quasi-particles in the ruthenate compound. These non- Abelian 
quasi-particles can be used as the states to be manipulated in a topological quantum 
computer. The biggest practical obstacle for quantum computation are the errors 
which may occur during a calculation due to materials being not ideal. Topological 
quantum computers minimize this source of error because they carry out operations 
by braiding the non- Abelian quasi-particles in a Hilbert-subspace containing degen- 
erate ground states. Due to an energy gap between this subspace and the rest of the 
Hilbert space it virtually decouples from all local perturbations |fT31 . Another con- 
firmed and well-studied condensed-matter example for an emerging p-wave struc- 
ture is superfluid 3 He-A JT6). 
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This section is a primer on the subject of spontaneous symmetry breaking, super- 
conductivity and superfluidity in the holographic context of the gauge/gravity corre- 
spondence. Little formalism is used, while we introduce all the necessary concepts. 



2.1 Basics of superconductivity & our field theory idea 

Let us review the essential concepts of superconductivity and understand how to 
build a p-wave superconductor. We are going to need this knowledge in order to ap- 
preciate the fact that our holographic setup reproduces this behavior in great detail. 



Table 1 Nomenclature of superconducting states 



Orbital anj 


mlar Name 


Parity of spatial Spin state 


momentum 




part 







s-wave 


even 


singlet 


1 


p-wave 


odd 


triplet 


2 


d-wave 


even 


singlet 



Superconductor basics & the p-wave Superconductivity is the phenomenon asso- 
ciated with infinite dc conductvity in materials at low temperatures. It is caused by 
the formation of a charged condensate in which directed currents do not experience 
resistivity. A defining criterion for superconductivity is the Meissner-Ochsenfeld 
effect described below. In conventional superconductors the superconducting con- 
densate consists of electron pairs called Cooper pairs. So there are two simultane- 
ous steps: the fermionic electrons have to pair up to form bosonic Cooper pairs, 
and these pairs do condense. This condensation happens in a second order phase 
transition (at vanishing magnetic field) when the temperature is lowered through its 
critical value T c . Due to the requirement for the fermionic state to be antisymmet- 
ric there are only certain symmetry combinations allowed for the two electron state 
describing a Cooper pair. As seen from table[T]the name "p-wave" superconductor 
refers to those pairs in which the relative orbital angular momentum between the 
two electrons is L = 1, the spatial part of the parity is odd and the spin state is a 
triplet. 

The mechanism pairing electrons in conventional superconductors is well de- 
scribed by a mean field theory approach and the microscopic BCS-theory . Recall 
that condensed matter systems are conveniently described in terms of lattices with 
many (about 10 23 ) sites. Then conventional BCS-theory tells us that the lattice is 
slightly deformed by the presence of an electron. This deformation can be de- 
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scribed by a quasi-particle excitation, a phonon. We can imagine the phonon to 
create a small potential well near the electron in which another electron can be 
caught. So lattice vibrations (phonons) mediate a weakly attractive interaction be- 
tween the electrons which then form pairs. Conventional BCS-theory with phonons 
is only valid at low temperatures because around 2QK the phononic lattice vibra- 
tions caused by the temperature already destroy the weakly attractive interaction 
between the electrons. Note that BCS-theory itself does not depend on the origin of 
the attractive interaction. 

A superconductor is simply a charged superfluid. The crucial defining property 
for any kind of superconductivity or superfluidity is that a symmetry is sponta- 
neously broken. In superfluids this symmetry is global while in superconductors it 
is local, i. e. a gauge symmetry. Therefore in superconductors there are a few addi- 
tional effects related to the gauge symmetry and the corresponding gauge field. But 
besides that those two phenomena are very similar. Especially in both cases there is 
a Goldstone boson created for each spontaneously broken symmetry of the describ- 
ing field theory. For a broken global symmetry this Goldstone boson survives and is 
visible in the spectrum as a hydrodynamic mode ifTTl . For a broken local symmetry 
however the Goldstone boson is eaten by the gauge field which couples to the charge 
belonging to the broken symmetry. In superconductors this causes the gauge boson, 
i. e. the photon for the broken electromagnetic t/(l) to become massive. Since these 
heavy photons can travel only an exponentially small distance, the electromagnetic 
interaction becomes short-ranged. Therefore magnetic fields, which can be thought 
of as consisting of photons, can only penetrate the system up to a certain distance, 
the penetration depth. This is called Meissner-Ochsenfeld effect and it is a defining 
criterion for superconductivity. In the Anderson-Higgs mechanism particles acquire 
a mass by the same mechanism. Thus it is sometimes described as the superfluidity 
of the vacuum. See lTT8l for a more precise review. 

There is a class of experimentally well-studied but theoretically less understood 
unconventional superconductors, such as copper or ruthenate compounds. Some of 
these materials show superconducting phases at high temperature^ up to 138K. 
The conventional BCS-theory does not apply to such high temperatures as men- 
tioned above. So the biggest mystery remains to understand the pairing mecha- 
nism of electrons in these high temperature superconductors. Pairing and condensa- 
tion need not occur at the same temperature here. Most important for the applica- 
tion of gauge/gravity duality: In these unconventional superconductors the coupling 
strength of the electrons to each other is strong. Thus these are experimentally ac- 
cessible systems governed by strongly coupled field theory. 

Another clear signature for superconductivity is of course an infinite dc conduc- 
tivity. Together with that there is a conductivity gap in the frequency-dependent 
conductivity. This is caused by the fact that the conductivity at small frequencies or 



4 These high temperature superconductors in general realize a d-wave structure. However, 
Sr2RuC>4 belongs to a distinct class of unconventional superconductors and is p-wave supercon- 
ducting at low temperatures around 2K. The pairing mechansim is not microscopically understood. 
The simplified argument is that the phonon-interaction of conventional Cooper pairing is isotropic, 
thus not providing an anisotropic p-wave structure. 
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energies vanishes until there is enough energy to break up one Cooper pair. From 
that energy on the material is a normal conductor with individual electrons being 
the charge carriers. In unconventional superconductors there is a surprising phe- 
nomenon called pseudo-gap. This means that the experiments carried out in un- 
conventional superconductors show a gap in the conductivity even at and above 
the transition temperature where the superconducting condensate forms. Inside this 
pseudo-gap the conductivity does not drop to zero but to a small finite value. 

How to build a field theory with p-wave superconductivity As stressed above 
the crucial thing to do in order to get a superconductor or superfluid is to sponta- 
neously break a symmetry. We are going to accomplish this by allowing our system 
to develop a charged condensate. 

Let us assume a particle physics point of view for a while. In these notes we will 
focus on a 3+1 - dimensional = 2 supersymmetric SU (N c ) Yang-Mills theory at 
temperature T, consisting of a jV = 4 gauge multiplet as well as Nf massive jY = 2 
supersymmetric hypermultiplets (y/j, 0;). The hypermultiplets give rise to the flavor 
degrees transforming in the fundamental representation of the gauge group. The 
action is written down explicitly for instance in |[T9l . In particular, we work in the 
large N c limit with Nf <C N c at strong coupling, i. e. with A >• 1 , where A = g\ M N c 
is the 't Hooft coupling constant. In the following we will consider only two flavors, 
i. e. Nf = 2. The flavor degrees of freedom are called u and d. If the masses of the 
two flavor degrees are degenerate, the theory has a global 1/(2) flavor symmetry, 
whose overall U(1)b subgroup can be identified with the baryon number. 

In the following we will consider the theory at finite isospin chemical potential 
jU, which is introduced as the source of the operator 

Jq<* VAC7 3 yoV+0CF 3 do0 =n u -n d , (1) 

where n u u is the charge density of the isospin fields, (0„, 0^) = and ( \j/ u , = \y. 
a 1 are the Pauli matrices. A non-zero vev (Jq) introduces an isospin density as 
discussed in l20l . The isospin chemical potential ji explicitly breaks the U{2) ~ 
U{\)b x SU(2)i flavor symmetry down to U(1)b x 17(1)$, where U( 1)3 is generated 

Fig. 1 Phase diagram of our 
field theory depending on the 
temperature parametrized by 
m <x T and on the isospin 
chemical potential fl 1 scaled 
by the quark mass M q . At low 
temperatures and chemical 
potentials a phase of stable 
mesons forms which melts at 
higher fl 1 and ]iT 1 . Above 
a critical density we find a 
superconducting or superfluid 
phase with a vector order 
parameter. 
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by the unbroken generator <7 3 of the SU(2)[. Under the Z7 (1)3 symmetry the fields 
with index u and d have positive and negative charge, respectively. 

However, the theory is unstable at large isospin chemical potential 11201 . The 
new phase is sketched in figure [T] We show in this lecture (see also Q), that the 
new phase is stabilized by a non-vanishing vacuum expectation value of the current 
component 

j\ °< y/a 1 73 y + a 1 d 3 <f> = f u 73 y/ c i + f d 7? y/„ + bosons . (2) 

This current component breaks both the 5(9(3) rotational symmetry as well as the 
remaining Abelian U(l)j, flavor symmetry spontaneously. The rotational SO (3) is 
broken down to 5(9(2)3, which is generated by rotations around the x 3 axis. Due 
to the non-vanishing vev for j\ , flavor charged vector mesons condense and form 
a superfluid. Let us emphasize that we do not describe a color superconductor on 
the field theory side, since the condensate is a gauge singlet. Figure Q] shows a "not 
accessable" parameter region in which we get divergent quantities. This is due to 
the fact that at such large charge densities we would need to take in account the 
backreaction of the D7-branes on the AdS geometry. 

In a condensed matter context our model can be considered as a holographic 
p-wave superconductor in the following way. The global t/(l) 3 in our model is 
the analog of the local i/(l) em symmetry of electromagnetic interactions. So far 
in all holographic models of superconductors the breaking of a global symmetry 
on the field theory side is considered. In our model, the current 7 3 corresponds 
to the electric current J em . The condensate breaks the J/(l) 3 spontaneously. 
Therefore it can be viewed as the superconducting condensate, which is analogous 
to the Cooper pairs. Since the condensate transforms as a vector under spatial 
rotations, it is a p-wave superconductor. - Strictly speaking, for a superconductor 
interpretation it would be necessary to charge the superfluid, i. e. gauge the global 
C/(l)3 symmetry which is broken spontaneously in our model. A spontaneously 
broken global symmetry corresponds to a superfluid. However, as we mentioned 
before many features of superconductivity do not depend on whether the t/(l) 3 is 
gauged. One exception to this is the Meissner-Ochsenfeld effect. To generate the 
currents expelling the magnetic field, the U(\)i symmetry has to be gauged. This 
matter is discussed further below in section 15.41 where we will be able to see the 
onset of the Meissner-Ochsenfeld effect within our holographic setup. 



2.2 Holographic realization 

Holographic superconductors have first been studied in ET1 l22l 171 l8l. The initial 
idea presented in i2T1 was that the Abelian Higgs model coupled to gravity with a 
negative cosmological constant provides a charged scalar condensate near but out- 
side a charged black hole horizon. The charged condensate spontaneously breaks 
the Abelian gauge symmetry of the theory. Later studies revealed that this break- 
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ing also occurs in setups with neutral black holes but with a negative mass for the 
charged scalar (7) . The basic idea for a holographic p-wave superconductor has been 
outlined in j8). Let us review the basic ideas. 

Holographic superconductor basics In order to spontaneously break a symmetry 
we need a gravity theory with a gauge symmetry. Note that this is not the usual 
gauge symmetry SU (N) of the correspondence where N — > °°, but an additional 
one with a finite rank, for example an Abelian 1/(1). Furthermore we need a black 
hole background in order to introduce finite temperature in the dual field theory, 
see J4] for details. Most importantly we need a charged condensate hovering in 
the bulk over the horizon. To be more precise, we need a bulk field charged 
under the gauge symmetry which is to be broken. <p has to have the typical ex- 
pansion = Normalizable + fynon-normalizable + ■■■ "ear the AdS boundary, but with 

0non-normaiizabie = 0. Why do we require this particular structure? Recall that in the 
gauge/gravity correspondence the normalizable mode is identified with the field the- 
ory expectation value of the operator dual to the gravity field <j>. This is our 
condensate in the field theory, which we want to be non-zero. The normalizable 
mode on the other hand is dual to a source in the field theory. Therefore we want 
it to vanish since it would break the symmetry explicitly. The role of the field <j> 
providing the condensate could be played for example by a charged scalar or by one 
component of a gauge field. These concepts are illustrated in the following example. 

Example: Consider the bottom-up toy model p-wave superconductor introduced 
in $8$. There we have an Einstein Yang-Mills theory with a negative cosmological 
constant 



with F being the field strength of an 51/(2) gauge field A. This SI/ (2) in this case 
is the gauge group that we want to break. Actually we break only aU(l) subgroup 
of it. But let us not worry about the details now. Only note that this theory is placed 
in an AdS^ charged black hole background. Our gauge field A now plays the role 
of the field (j) providing the condensate. The operator dual to the gauge field A is 
the electromagnetic current = J^. The boundary behavior of the gravity field's 
components according to the equations of motion derived from equation (0 is given 



with the radial AdS coordinate r. In principle we could allow more non-vanishing 
components but this combination turns out to be both sufficient and consistent. As 
usual in thermal AdS/CFT the temporal component A t introduces a chemical poten- 
tial ll in the dual field theory. This chemical potential sources the corresponding 
operator J' = d which is simply the charge density of the SU(2) charge in the field 




(3) 



by 




(4) 
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theory at the boundary. This explicitly breaks SU (2) — > U(l)^. However the spatial 
component A\ spontaneously breaks the remaining t/(l)3 with the condensate (J^). 

How to build a gravity dual to p-wave superconductivity As mentioned several 
times before, we need a vector condensate for our p-wave superconductor. Con- 
veniently the previous example already showed us what the structure for the dual 
gravity theory has to be in order to realize a vector condensate. Unfortunately that 
example has not been derived from string theory. Probably the most difficult task in 
building a holographic superconductor is finding a gravity setup with all the necces- 
sary features, which is actually stable and thermodynamically favored. We are going 
to see that the system of intersecting D3 and D7 branes provides exactly such a sta- 
ble configuration. Since we are going to review Dp/Dq- brane systems below, let 
us here spot only those features which are of importance to superconductivity. The 
SU (2) "flavor" gauge group which we are going to break (partly spontaneously) is 
created by using two coincident D7 branes. By introducing a gauge field A^ living 
on these D7-branes and giving its temporal component a non-trivial profile in one of 
the flavor directions, we introduce a chemical potential in the dual field theory. This 
is completely analogous to introducing the non-trivial A^ in the previous example. 
All this is going to take place in the background of a stack of N c black D3-branes. 
They introduce the temperature into the dual field theory. Finally, the gravity field 
which will break the residual flavor symmetry is going to be a spatial component 
of the gauge field just as Aj. in the previous example. Our analysis of the thermo- 
dynamic potentials is going to show that this phase is thermodynamically preferred 
over the phase without the symmetry-breaking condensate. Furthermore it is stable 
against all obvious gauge field fluctuations. 



3 Holographic setup 

In this section we carry out exactly the D3-D7 brane construction outlined in section 
12.21 The result will be a gravity setup being holographically dual to a p-wave super- 
fluid/superconductor. But let us first review how to add flavor to the gauge/gravity 
correspondence. 



3.1 Flavor from intersecting branes 

Let us imagine for this subsection that we want to use gauge/gravity in order 
to model QCD or the quark gluon plasma state of matter produced at the RHIC 
Brookhaven heavy-ion collider. The original AdS/CFT conjecture does not include 
matter in the fundamental representation of the gauge group but only adjoint matter. 
In order to come closer to a QCD-like behavior one can therefore investigate how 
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to incorporate quarks and their bound states in this section. We focus on the main 
results of 11231 and 11241 . however for a concise review the reader is referred to |fl9l . 

Since AdS/CFT has been discovered a lot of modifications of the original con- 
jecture have been proposed and analyzed. This is always achieved by modifying the 
gravity theory in an appropriate way. For example the metric on which the gravity 
theory is defined may be changed to produce chiral symmetry breaking in the dual 
gauge theory l25l l26l . Other modifications put the gauge theory at finite temper- 
ature and produce confinement IT271 . Besides the introduction of finite temperature 
the inclusion of fundamental matter, i.e. quarks, is the most relevant extension for us 
since we are aiming at a qualitative description of strongly coupled QCD effects at 
finite temperature. These effects are similar to the ones observed at the RHIC heavy 
ion collider. 

Adding flavor to AdS/CFT The change we have to make on the gravity side in 
order to produce fundamental matter on the gauge theory side is the introduction 
of a small number Nf of D7-branes. These are also called probe branes since their 
backreaction on the geometry originally produced by the stack of N D3-branes is 
neglected. Strings within this D3/D7-setup now have the choice of starting (ending) 
on the D3- or alternatively on the D7-brane. Note that the two types of branes share 
the four Minkowski directions 0,1,2,3 in which also the dual gauge theory will 
extend on the boundary of AdS as visualized in figure [2] The configuration of one 
string ending on N coincident D3-branes produces an SU (N) gauge symmetry of 
rotations in color space. Similarly the Nf D7-branes generate a U(Nf) flavor gauge 
symmetry. We will call the strings starting on the stack of Dp-branes and ending 
on the stack of Dg-branes p — q strings. The original 3 — 3 strings are unchanged 
while the 3 — 7- or equivalently 7 — 3 strings are interpreted as quarks on the gauge 
theory side of the correspondence. This can be understood by looking at the 3 — 3 
strings again. They come in the adjoint representation of the gauge group which can 
be interpreted as the decomposition of a bifundamental representation (A^ 2 — 1)0 
1 = N <S)N. So the two string ends on the D3-brane are interpreted as one giving 
the fundamental, the other giving the anti-fundamental representation in the gauge 
theory. In contrast to this the 3 — 7 string has only one end on the D3-brane stack 
corresponding to a single fundamental representation which we interpret as a single 
quark in the gauge theory. 

We can also give mass to these quarks by separating the stack of D3-branes from 
the D7-branes in a direction orthogonal to both branes. Now 3 — 7 strings are forced 
to have a finite length L which is the minimum distance between the two brane 
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Fig. 2 Coordinate directions in which the Dp-branes extend are marked by 'x'. D3- and D74}ranes 
always share the four Minkowski directions and may be separated in the 8,9-directions which are 
orthogonal to both brane types. 
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stacks. On the other hand a string is an object with tension and if it assumes a mini- 
mum length, it needs to have a minimum energy being the product of its length and 
tension. The dual gauge theory object is the quark and it now also has a minimum 
energy which we interpret as its mass M q = L/(2na'). 

The 7 — 7 strings decouple from the rest of the theory since their effective cou- 
pling is suppressed by Nf/N. In the dual gauge theory this limit corresponds to 
neglecting quark loops which is often called the quenched approximation. Never- 
theless, they are important for the description of mesons as we will see below. 

Let us be a bit more precise about the fundamental matter introduced by 3 — 7 
strings. The gauge theory introduced by these strings (in addition to the original 
setup) gives a jY = 2 supersymmetric U(N) gauge theory containing Nf funda- 
mental hypermultiplets. 

D7 embeddings & meson excitations Mesons correspond to fluctuations of the 
D7-brane^f| embedded in the AdS$ x S 5 -background generated by the D3-branes. 
From the string-point of view these fluctuations are fluctuations of the hypersurface 
on which the 7 — 7 strings can end, hence these are small oscillations of the 7 — 7 
string ends. The 7 — 7 strings again lie in the adjoint representation of the flavor 
gauge group for the same reason which we employed above to argue that 3 — 3 
strings are in the adjoint of the (color) gauge group. Mesons are the natural objects 
in the adjoint flavor representation. Vector mesons correspond to fluctuations of the 
gauge field on the D7-branes. 

Before we can examine mesons as D7-fluctuations we need to find out how the 
D7-branes are embedded into the 10-dimensional geometry without any fluctua- 
tions. Such a stable configuration needs to minimize the effective action. The ef- 
fective action to consider is the world volume action of the D7-branes which is 
composed of a Dirac-Born-Infeld and a topological Chern-Simons part 

Sdi = -T D1 ,/d 8 cj e -* yj- det {P[g + B} al3 + (2na')F ap } (5) 

+ { ^^-T D7 JP[C 4 ]AFAF. (6) 

The preferred coordinates to examine the fluctuations of the D7 are obtained from 
the standard AdS coordinates 

ds 2 = ^dp 2 + |J(-dr + dx 2 ), (7) 

with the AdS radius R and the dimensionful radial AdS coordinate p. 
Exercise Show that the standard AdS metric (0 transforms to 

r 2 R 2 

ds 2 = -^dx 2 + — (dp 2 + p 2 di2 3 2 + dw 5 2 + dw 6 2 ) , (8) 



5 To be precise the fluctuations correspond to the mesons with spins 0, 1/2 and 1 I24II28I . 
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under the transformation p = w\ + . . . + W4 , r = p + w$ + W(, , where x is a 
four vector in Minkowski directions 0, 1,2,3 and R is the AdS radius. The coordi- 
nate r is the radial AdS coordinate while p is the radial coordinate on the coincident 
D7 -branes. 

Let us follow l24l : For a static D7 embedding with vanishing field strength F on 
the D7 world volume the equations of motion are 



(9) 



where wsfi denotes that these are two equations for the two possible directions of 
fluctuation. Since (O is the equation of motion of a supergravity field in the bulk, 
the solution near the AdS boundary takes the standard form with a non-normalizable 
and a normalizable mode, or source and expectation value respectively 

c 

w 56 = m+—f + ... , (10) 
P 2 

with m being the quark mass acting as a source and c being the expectation value 
of the operator which is dual to the field W5 6- While c can be related to the scaled 
quark condensate c « (qq)(2na') 3 . 

If we now separate the D7-branes from the stack of D3-branes the quarks become 
massive and the radius of the 5 3 on which the D7 is wrapped becomes a function of 
the radial AdS coordinate r. The separation of stacks by a distance m modifies the 
metric induced on the D7 P[g] such that it contains the term R 2 p 2 /(p 2 + m 2 )dQ^ 2 . 
This expression vanishes at a radius p 2 = r 2 — m 2 = such that the S 3 shrinks to 
zero size at a finite AdS radius. 

Fluctuations about these w$ and w& embeddings give scalar and pseudoscalar 
mesons. We take 

h>5 = + 2na'x , W6=m+2na'(p (11) 

After plugging these into the effective action (0 and expanding to quadratic order 
in fluctuations we can derive the equations of motion for <p and %. As an example 
we consider scalar fluctuations using an Ansatz 

(p = $(p)e ikx &,(S 3 ), (12) 

where ^(S 3 ) are the scalar spherical harmonics on the S 3 , solves the radial part 
of the equation and the exponential represents propagating waves with real momen- 
tum k. We additionally have to assume that the mass-shell condition 

M 2 = k 2 (13) 

is valid. Solving the radial part of the equation we get the hypergeometric func- 

2 

tion <j) oc F(—a, — 0C + 1+ 1; (1 + 2); z ^ r ) and the parameter 
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1- Jl-k 2 R 4 /m 2 
a = ^ ( i4) 

summarizes a factor appearing in the equation of motion. In general this hypergeo- 
metric function may diverge if we take p — > °°. But since this is not compatible with 
our linearization of the equation of motion in small fluctuations, we further demand 
normalizability of the solution. This restricts the sum of parameters appearing in the 
hypergeometric function to take the integer values 

n=a-l-l, n = 0,1,2,.... (15) 

With this quantization condition we determine the scalar meson mass spectrum to 
be 

2m r- — - 

M(pseudo) S calar = -j V (« + ' + 1 ) (« + ' + 2 ) , (16) 

where « is the radial excitation number found for the hypergeometric function. Sim- 
ilarly we can determine pseudoscalar masses and find the same formula[T6l For vec- 
tor meson masses we need to consider fluctuations of the gauge field A appearing in 
the field strength F in equation ((5]). The formula for vector mesons (corresponding 
to e.g. the p -meson of QCD) is 

2m 



M v = -jpy/(n + l + l)(n + l + 2). (17) 

Note that the scalar, pseudoscalar and vector mesons computed within this frame- 
work show identical mass spectra. Further fluctuations corresponding to other 
mesonic excitations can be found in Il24ll28l . 

Brane embeddings at finite temperature In order to get a finite temperature in 
the dual field theory, the gravity theory needs to be put into a black hole or black 
brane background geometry. It was found in l26l |29l [30l that at finite temperature 
our Nf probe flavor branes can be embedded in two distinct ways. There are high 
temperature configurations called black hole embeddings in which part of the brane 
falls into the black hole horizon. On the other hand there are low-temperature con- 
figurations called Minkowski embeddings in which the brane stays outside the black 
hole horizon. These two configurations are separated by a geometric transition, i. e. 
the configuration in which the brane just barely touches the black hole horizon. This 
geometric transition corresponds to the meson melting transition for the fundamen- 
tal matter of the dual field theory. Note that the adjoint matter of this field theory is 
always deconfined in this setup. 

At finite charge density however there is only one kind of embedding and that 
is the black hole embedding. The heuristic argument is that introducing a finite 
charge on the brane there have to be field lines for the associated field strength. 
These lines have to end somewhere. Our setup has rotational symmetry in the spatial 
directions. Imagine the radial AdS coordinate and field lines running along it. If they 
are supposed to end somewhere, there has to be a horizon. Otherwise they would all 



Flavor Superconductivity & Superfluidity 



15 



meet in the origin at p = 0. Note that field lines ending at the horizon means, that we 
can interpret the horizon as being charged. In these lecture notes we will exclusively 
deal with non-vanishing density and thus only encounter black hole embeddings. 



3.2 Background and brane configuration 

We aim to have a field theory dual at finite temperature. This is holographically ac- 
complished by placing the gravity theory in a black hole or black brane background. 
Here the black hole's Hawking temperature can be identified with the field theory 
temperature. We consider asymptotically AdSs x S 5 space-time. The AdSs x S 5 ge- 
ometry is holographically dual to the jV = 4 Super Yang-Mills theory with gauge 
group SU (N c ). The dual description of a finite temperature field theory is an AdS 
black hole. We use the coordinates of OP to write the AdS black hole background 
in Minkowski signature as 

d * 2 = ^ (-£dt 2 +fdA + (^X(dp 2 + p 2 dQ 2 ), (18) 



2R 2 \ f J \P 

with dQ,^ the metric of the unit 5-sphere and 

/(p)=i-4, /(p)=i+4, d9) 

r r 

where R is the AdS radius, with R 4 = 4-ng s N c a' 2 = 2X a' 2 . 

Exercise: Show that the metric M8$ can be obtained from the standard finite tem- 
perature AdS metric 

(kTR^ 2 R 2 
ds 2 = y J [-/( M )df 2 + dx 2 ] + 4 ^^^ dM 2 +^ 2 dX2 5 2 ; /(m) = 1-m 2 , (20) 



by the transformation (rop) 2 = 2r 2 + \J r 4 — ro 4 , with r being the original radial 
AdS coordinate and ro the location of the black hole horizon. 

Exercise: The temperature of the black hole given by M8t may be determined by 
demanding regularity of the Euclidean section. Show that it is given by 

r =H- (21) 

In the following we may use the dimensionless coordinate p = p / p#, which 
covers the range from the event horizon at p = 1 to the boundary of the AdS space at 
p — > °°. To include fundamental matter, we embed Nf coinciding D7-branes into the 
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ten-dimensional space-time as illustrated in figure [2] These D7-branes host flavor 
gauge fields with gauge group U (Nf). To write down the DBI action for the D7- 
branes, we introduce spherical coordinates {r,Qi} in the 4567-directions and polar 
coordinates {L, 0} in the 89-directions [31]. The angle between these two spaces is 
denoted by 9 (0 < G < n/2). The six-dimensional space in the 456789-directions is 
given by 

dp 2 + p 2 dX2 5 2 = dr 2 + r 2 di2 3 2 + dL 2 + L 2 d<j) 2 = dp 2 + p 2 (d 2 + cos 2 0d<j) 2 + sin 2 6dQ%) , 

(22) 

where r = p sin 0, p 2 = r 2 +L 2 and L = p cos 0. D3- and D7-branes always share the 
four Minkowski directions and may be separated in the 8,9-directions which are or- 
thogonal to both brane types. That separation is dual to the mass of the fundamental 
flavor fields in the dual gauge theory, i. e. the quarks. 

Due to the 5(9(4) rotational symmetry in the 4567 directions, the embedding of 
the D7-branes only depends on the radial coordinate p. Defining % = cos0, we 
parametrize the embedding by % = x(P) an d choose 0=0 using the 5(9(2) sym- 
metry in the 89-direction. The induced metric G on the D7-brane probes is then 

p 2 l_ x 2 

(23) 

The square root of the determinant of G is given by 

= ^V//(l - r) V 1 - X 2 + P 2 ( d pX) 2 , (24) 

where is the determinant of the 3-sphere metric. 

As in l20l we introduce a SU (2) isospin chemical potential jj. by a non-vanishing 
time component of the non-Abelian background field on the D7-brane. The genera- 
tors of the SU(2) gauge group are given by the Pauli matrices a'. Due to the gauge 
symmetry, we may rotate the flavor coordinates until the chemical potential lies in 
the third flavor direction, 

H= lim A 3 (p). (25) 

This non-zero gauge field breaks the SU(2) gauge symmetry down to £7(1)3 gen- 
erated by the third Pauli matrix a 3 . The spacetime symmetry on the boundary is 
still 5(9(3). Notice that the Lorentz group 5(9(3, 1) is already broken down to 5(9(3) 
by the finite temperature. In addition, we consider a further non-vanishing back- 



^ 2(G) = 2^("7 df2+/dx2 ) 



coord, names 




/i, v, . . . 

Fig. 3 This figure summa- indices i j 

rizes our coordinates and 12 

indices. 
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ground gauge field which stabilizes the system for large chemical potentials. Due 
to the symmetry of our setup we may choose A^dx 3 ^ 1 to be non-zero. To obtain an 
isotropic configuration in the field theory, this new gauge field A3 only depends on 
p . Due to this two non-vanishing gauge fields, the field strength tensor on the branes 
has the following non-zero components, 



^3 = 




= d p A\ 


(7 — 7strings) , 


(26) 






= dpAl 


(3 — 7strings) . 


(27) 


F 03 = - 


" F 30 = 


- -J-A 3 A 1 

' VI 03 


(interaction) . 


(28) 



The labels behind those equations refer to the sort of strings which generate the 
corresponding gauge fields. The field strength Fq 3 can be understood as an interac- 
tion term between 7-7 and 3-7 strings. We derive this interpretation in section loTTI 



3.3 DBI action and equations of motion 

In this section we calculate the equations of motion which determine the profile 
of the D7-brane probes and of the gauge fields on these branes. A discussion of 
the gauge field profiles, which we use to give a geometrical interpretation of the 
stabilization of the system and the pairing mechanism, may be found in section loTTl 
The DBI action determines the shape of the brane embeddings, i. e. the scalar 
fields (j), as well as the configuration of the gauge fields A on these branes. We 
consider the case of Nf = 2 coincident D7-branes for which the non-Abelian DBI 
action reads 



S DB i = -7b7Str /d 8 £ Vdetg 
with 



dalP^E^+E^Q- 1 -8) ij E jv ] +2xa'F ab 



(29) 



Q'j = 8'j + i2Ka'[4>\ <P k ]E kj (30) 

and P a h the pullback to the Dp-brane, where for a D/j-brane in d dimensions we have 
H, v =0,...,(d- 1), a,b = 0,...,p, i, j = (p + 1),. . . , (d- 1), E MV = guv +B^ V . 
In our case we set p = 7, d = 10, B = 0. As in ll20l we can simplify this action 
significantly by using the spatial and gauge symmetries present in our setup. The 
action becomes 

Sdbi = -Ten J d 8 £ Str A /|det(G + 2ffa'F)| (31) 
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1+G (,0 G 44 (Fi (a 3 ) +G 33 G [Fp 3 \ (a 1 ) 



(32) 



where in the second line the determinant is calculated. Due to the symmetric trace, 
all commutators between the matrices a' vanish. It is known that the symmetrized 
trace prescription in the DBI action is only valid up to fourth order in a' IT331 [34ll . 



However the corrections to the higher order terms are suppressed by N7 1 1331 (see 
also |36|). Here we tise two different approaches to evaluate the non-Abelian DBI 
action OTb . First, we modify the symmetrized trace prescription by omitting the 
commutators of the generators a' and then setting (cr 1 ) 2 = 1 (see subsection 13. 3. II 
below). This prescription makes the calculation of the full DBI action feasible. This 
prescription is not verified in general but we obtain physically reasonable results as 
discussed in section I5TI and [5721 Second, we expand the non-Abelian DBI action to 
fourth order in the field strength F (see subsection l3.3.21 i. Here it should be noted 
that in general the higher terms of this expansion need not be smaller than the lead- 
ing ones. However, we again get physical results in our specific case which confirm 
this approach. We further motivate the validity of our two approaches below. 



3.3.1 Adapted symmetrized trace prescription 

Using the adapted symmetrized trace prescription defined above, the action becomes 



Sdbi = -TwNf d^V^G 



1 (^ ) 2 +g 33 g 44 (f; 3 ) 2 +g 00 g 33 (4) 2 



T nnN f , ^ p3//(1 _^ )r(p ^ A) 



(33) 



with 



1 - X 2 + P 2 (d P X) 2 - 2 ^(l-X 2 ) (d p Al) 2 + j (1 - X 2 ) (d p A\] 

-j^(i-x 2 +P 2 (d P x) 2 )(m) 



(34) 



where the dimensionless quantities p = p/pn and A = (2na')A/ pn are used. To 
obtain first order equations of motion for the gauge fields which are easier to solve 
numerically, we perform a Legendre transformation. Similarly to ||3"T1 |20l we cal- 
culate the electric displacement and the magnetizing field p\ which are given by 
the conjugate momenta of the gauge fields Al and Al, 
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In contrast to BTl [37l l38l l20l . the conjugate momenta are not constant any more 
but depend on the radial coordinate p due to the non-Abelian term AqA^ in the DBI 
action. For the dimensionless momenta /5q and p\ defined as 

(36) 



Ina'NfTo-jpfj 
we get 

Po ~ 2fY( P>X) A) ' P3 ~ 2T(p^,A)) • (3?) 
Finally, the Legendre-transformed action is given by 



'DBI 



>DB1 



(d p A 



3\ SS D m 



(Pi a 1 \ Ss ° B1 



(38) 



Id^ P 3 ff(l-x 2 )V^-X 2 + P 2 (d pX ) 2 V(p,X,lp), (39) 



with 



V(p,X,A,/>) = 



1 27 (A?A 1 ) 2 



X 1 + 



s n 2 



p 6 P(i-x 2 ) 3 p 6 // 2 (i-x 2 ) 3 



(40) 



Then the first order equations of motion for the gauge fields and their conjugate 
momenta are 



d P A o 



d p A\ 



73 _ 2f v /l-x 2 +p 2 (dp X ) 2 - 3 



<5pP0 = 



p 3/2 (1 _ x 2)2 P0 W (P'^= A 'P), 

/(l-X 2 )Vl-X 2 +P 2 (^X) 2 ^ 2 ,'11x2 13 
2^p/W(p, Z ,A,;5) 



.1 _ f(l-X 2 )^l-X 2 +P 2 (d P x) 2 c 2 ,7 3 \2 ri 



with 



<5p/ 5 3 



A o) ^3 ' 



(41) 
(42) 
(43) 
(44) 

(45) 



\ 1t pW7F p 6 // 2 (i-z 2 ) 3 

For the embedding function % we get the second order equation of motion 
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p 5 ff(l-X 2 )(d P X) 
^/l-X 2 + P 2 (d P X) 2 



V 



p 3 ffx 



y/l-X 2 +p 2 (dpX) 2 



[3(l- X 2 )+2p 2 (dpX) 2 ]V 



2A(i-x 2 +pHd pX f) r, 3 , 2 p , 



(46) 



We solve the equations numerically and determine the solution by integrating the 
equations of motion from the horizon at p = 1 to the boundary at p = °o. The initial 
conditions may be determined by the asymptotic expansion of the gravity fields near 
the horizon 



4 3 
A 



^3 



Co 



Po = c o 

Pi = 



X=Xo 



y/(i-z&)>+4 

b +<?((p_l)3) , 



(p-l) 2 +^((p-l) 3 ) 



+ ^((p-l) 3 ), 
3Zo(l-Zo) 3 



m-xlY+4 



0{{p-l? 



.(p-l) 2 +^((p-l) 3 ) 



(47) 

(48) 

(49) 
(50) 

(51) 



where the terms in the expansions are arranged according to their order in p — 1 . For 
the numerical calculation we consider the terms up to sixth order in p — 1 . The three 
independent parameter bo, co and Xo ma Y be determined by field theory quantities 
defined via the asymptotic expansion of the gravity fields near the boundary, 



(52) 

(53) 
(54) 
(55) 
(56) 



According to the AdS/CFT dictionary, jj, is the isospin chemical potential. The pa- 
rameters d are related to the vev of the flavor currents J by 



A 3 o 


= P 








A\ 




A+° 
p~ 


(P 




Pi 


= 4 




• 




P\ 


= -d\ 






Hp- 4 


4n 2 p 2 


+ t 




X 




m c 

P P i 


-6 


(p- 4 ) 



2H4) 



d\ 



N f N c VXT 3 ' N f N c \/XT 3 
and m and c to the bare quark mass M q and the quark condensate (\jf\jf), 



(57) 
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2M a 8(w) 
™=^J-, c = - - , (58) 

respectively. There are two independent physical parameters, e. g. m and jx, in the 
grand canonical ensemble. From the boundary asymptotics d52l ). we also obtain that 
there is no source term for the current 7 3 ' . Therefore as a non-trivial result we find 
that the t/ (1)3 symmetry is always broken spontaneously. In contrast, in the re- 
lated works on p-wave superconductors in 2 + 1 dimensions l8l [39l , the spontaneous 
breaking of the 1/(1)3 symmetry has to be put in by hand by setting the source term 
for the corresponding operator to zero. With the constraint A3 |p^oo = and the two 
independent physical parameters, we may fix the three independent parameters of 
the near-horizon asymptotics and obtain a solution to the equations of motion. 



3.3.2 Expansion of the DBI action 



We now outline the second approach which we use. Expanding the action PTT ) to 
fourth order in the field strength F yields 



Sdbi = -TmNf / d 8 £ V^G 



2 8 



(59) 



where ^ consists of the terms with order i in F. To calculate the we use the 
following results for the symmetrized traces 



2a : Str 
4cj: Str 



= N f , 

= N f , Str[(c7')V) 2 



N f 

T' 



(60) 
(61) 



where the indices i,j are distinct. Notice that the symmetric trace of terms with un- 
paired <j matrices vanish, e. g. Str[<7'<7 7 ] =Nfd' j . The ^ are given in the appendix 

of El. 

To perform the Legendre transformation of the above action, we determine the 
conjugate momenta as in d35l l. However, we cannot easily solve these equations for 
the derivative of the gauge fields since we obtain two coupled equations of third 
degree. Thus we directly calculate the equations of motion for the gauge fields on 
the D7-branes. The equations are given in the appendix of (6). 

To solve these equations, we use the same strategy as in the adapted symmetrized 
trace prescription discussed above. We integrate the equations of motion from the 
horizon at p = 1 to the boundary at p = °° numerically. The initial conditions may 
be determined by the asymptotic behavior of the gravity fields near the horizon 



Al = a 2 {p-\) 2 + ff{(p-\) 
A l 3 = b o + 0((p-l) 3 ) , 



(62) 
(63) 
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*= * + 3 ffi^+> - 1)2 + "&-V 3 )- (64) 

For the numerical calculation we use the asymptotic expansion up to sixth order. 
As in the adapted symmetrized trace prescription, there are again three independent 
parameters a2,bo,Xo- Since we have not performed a Legendre transformation, we 
trade the independent parameter cq in the asymptotics of the conjugate momenta 
in the symmetrized trace prescription with the independent parameter «2 (cf. asymp- 
totics in equation (|47ii). However, the three independent parameters may again be 
determined in field theory quantities which are defined by the asymptotics of the 
gravity fields near the boundary 

<P 

Al = ii--± + 0{p 4 ) , (65) 

r 

d l 

A 3= -7^ + (P 4 ) > (66) 

m c , a-* 
X= p+p5 + *{P 4 )- ( 67 > 

The independent parameters ^i,d^,d\,m,c are given by field theory quantities as 
presented in ( fSTb and d58b . Again we find that there is no source term for the cur- 
rent /j, which implies spontaneous symmetry breaking. Therefore the independent 
parameters in both prescriptions are the same and we can use the same strategy to 
solve the equations of motion as described below i 



Exercise: In H401I only the leading order of the action ( 1591 ) quadratic in field fluc- 
tuations was considered. For a specific chemical potential an analytic solution with 
non-zero Aq and A3 can be found. Show that the analytic solution found there 
(adapted to ourAdS$ case) 

Al = fi(l-p 4 ), A\ = e {l f p4)2 , (68) 

indeed solves the equations of motion derived from the action expanded to quadratic 
order for the chemical potential 11 — 4. £ is a constant formed from the coupling 
constant and the vacuum expectation value of the dual current {J\). Note: This 
exercise requires some work. See also BJ^ for an application of this solution. 



4 D-brane thermodynamics & spectrum 



In this section we briefly review the results obtained for the thermodynamics and 
the spectrum of our setup with a baryonic, and later an isospin chemical potential. 
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Figure |4] shows the phase diagram of the field theory dual to the gravity setup we 
have constructed in the previous section [3] In this field theory we have introduced 
a baryonic chemical potential jj. q which is shown on the vertical axis scaled by the 
quark mass M q defined by equation d58l l. It is introduced by a non trivial back- 
ground gauge field solving the equation of motion derived from the DBI-action and 
asymptoting to the chemical potential near the boundary as lim p ^ PB A t (p) = ji q . 
The horizontal axis shows the scaled temperature T/M = m _1 . Recall that m is the 
asymptotic value of the D7-brane embedding near the AdS boundary. In other words 
it is the source term for a quark condensate, or the non-normalizable mode of the 
brane embedding function. Here however we use it merely as a temperature scale. 
The lines in the diagram in figure |4] are lines of equal baryon density d. At low 
temperature and chemical potential there is a triangle-shaped phase of zero density. 
Its diagonal borderline to the white region with finite baryon density is the location 
of the so-called meson-melting transition. This is the transition where the funda- 
mental matter melts, i. e. quark bound states, the mesons, acquire an increasing 
decay width becoming quasi-particles. In the grey phase we therefore have stable 
mesons with zero decay width. The deeper we go into the white phase away from 
the transition line, the mesons melt. We will also see this in the spectral functions 
computed below. As explained at the end of section 13.11 we are going to stay en- 
tirely in the white phase at finite temperature where the brane embeddings are of 
black hole type, which means that part of them falls into the black hole horizon (see 
section IXTT i. 

4.1.1 Correlator recipe 

We need to examine fluctuations of the fields in our setup in order to determine its 
spectrum and stability. There might be a fluctuation which is tachyonic and therefore 
could destabilize the whole system. For this purpose we quickly review how to 
compute real-time correlation functions in gauge/gravity. 



Fig. 4 The phase diagram in 
the canonical ensemble plot- 
ted against the variables of 
the grandcanonical ensemble. 
On the axes the scaled chem- 
ical potential \l q /M q , with 
the quark mass M q is shown 
versus the scaled tempera- 
ture T/M = mT 1 from (42). 
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Let us work along the example of a gauge field fluctuation a p. This appears in the 
action, in our case the Dirac-Born-Infeld action, in the field strength F = da + a A a. 
Note that a background gauge field A would appear in this field strength too, as we 
will see in later sections. For simplicity we consider Abelian gauge field fluctuations 
without background A, i. e. Fp V = dpa v — d v ap. The action then reads something 
like 



S~Jd*Zy/P\g llv ]+F flv , (69) 

with the P being the pullback of the metric g to the flavor brane. Expanding this 
action to quadratic order in fluctuations a, we get the linearized equation of motion 
for example for the spatial fluctuation a y which in Fourier space looks like this 

2 , d p [y/=g^gPP] . , g" 2 



= dp a y H — -dpdy + —w a y , (70) 

v~ 88 ■ ■ " 8 

where w is the dimensionless frequency of the fluctuation. Also for other fields we 
end up with a second order differential equation that we need to solve. Usually these 
equations have singular (at the horizon p#) coefficients which need to be regular- 
ized by an appropriate ansatz. In order to find the most singular behavior solving 
this equation, we plug the ansatz (p — Ph)^ into equation ( TTOb and expand. The 
leading order is a quadratic equation which can be solved for j3 giving /3; n or j3 0Ut . 
Only one of those two solutions, j8; n describes a fluctuation falling into the horizon, 
the other one is outgoing. We discard the outgoing one because nothing is supposed 
to leave a classical black hole. j3 is sometimes called the indicial exponent. Now 

we plug a = (p - pt{)P m F(p), with F(p) = £ f„(p - pn)" being a regular func- 

tion of p, into the equation of motion d70l ). Note that there might be logarithmic 
terms present as explained in general for example in 1431 and discussed in detail in 
J4). Picking the ingoing wave fixes one of the two boundary conditions. The other 
boundary condition may be fixed by choosing the normalization of F(p), i. e. the 
value /o. All the higher /„ depend recursively on /o and the indicial exponent p\ 
Now we solve the equation for F analytically or numerically. The correlator is then 
obtained from the quadratic part of the on-shell action, which in our case has the 
structure S on -sheii ~ / dpB(p)ad p a, where B is a function of p depending on metric 
coefficients. The recipe developed in P4l |451 tells us to strip the boundary values 
a bdy Q ff f rom j-jjg gelds a = a bdy &/ and identify what remains of the integrand with 
the Green function at the boundary p = pg 

= lim B{p)^pd p ^ v . (71) 



Now we only need to plug in our solutions a = a My £/. More details on the analytic 
and numerical procedures are explained in |4j. 
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The spectral function ffl is obtained from the imaginary part of the Green function 
2% = — 2lmG R . It encodes the spectrum of our thermal field theory. In particular in 
figure|5] we are able to identify pretty stable quasi-particle excitations at a low finite 
temperature parametrize by Xo an d at finite baryon density d = 0.25. Lowering the 
temperature these quasiparticles approach the line spectrum ( fT71 ) indicated here as 
dashed vertical lines. The reason is that at lower temperature our theory restores 
its original supersymmetry. This also tells us that the vector quasi-particles we see 
in the thermal spectral function are vector mesons. Our mesons melt in the finite 
temperature, finite density phase as mentioned above in the discussion of the phase 
diagram H 

Quasi normal modes (QNM) The difference between the zero temperature line 
spectrum and the finite temperature spectrum of quasi-particle excitations lies in 
the nature of the corresponding eigenmodes. In the zero temperature case there is 
no black hole, thus no dissipation on the gravity side. The system has well-defined 
normal modes at real frequencies w £ R. At finite temperature however the (quasi) 
eigenfrequencies are complex w € C owing to the dissipation into the field theory 
plasma or in the dual gravity picture: dissipation into the black hole. The modes trav- 
eling with these (quasi) eigenfrequencies are called quasi-normal modes (QNM). 
We can roughly think of quasi-normal modes as being those solutions to the gravity 
fluctuation equations which vanish at the AdS boundary. These QNMs have been 
found to be identical to the poles in the field theory Green function. Therefore the 
location of the QNMs is closely related to the location of the peaks in the spectral 
function. At least some of the quasi-normal modes create quasi-particle peaks in 
the spectral function. In the zero temperature limit we can think of those QNMs 
as approaching the real frequency axis and reaching it in the limit, becoming real- 
valued. The corresponding quasiparticles become stable which means line-shaped 
in the spectrum. A more complete picture of QNMs is given in [46|. Quasi normal 
modes of our particular D3/D7 system are discussed in l47ll48l . 



Fig. 5 The spectral function 
ffi (in units of NfN c T 2 /4) at 
finite baryon density d versus 
dimensionless frequency w. 
At large Xo here the peaks 
approach the dashed drawn 
line spectrum given by G3. 



E 
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4.2 Isospin chemical potential 

We can equally well introduce a chemical potential with SU (2) -let's call it isospin- 
structure, represented by the Pauli matrices a'. Choosing the chemical potential 
to point along the third flavor direction this background boils down to having two 
copies of the Abelian background gauge field A t (p ) explored above in the following 

way 

However we get some interesting new signatures through the SU (2) structure. For 
example the spectrum shown in figure [6] shows a triplet splitting for our mesons. 
In particular we observe a splitting of the line expected at the lowest meson mass 
at w = 4.5360 (« = 0). The resonance is shifted to lower frequencies for £%xy and to 
higher ones for Myx, while it remains in place for & E i E i. The second meson res- 
onance peak (n = 1) shows a similar behavior. So the different flavor combinations 
propagate differently and have distinct quasi-particle resonances. This behavior is 
analogous to that of QCD's p meson which is vector meson and a triplet under the 
isospin SU (2) of QCD. Thus we have modeled the melting process of vector mesons 
in a quark gluon plasma at finite isospin density. 



4.3 Instabilities & the new phase 

What does all this have to do with our condensed matter motivation? The cru- 
cial thing is that the isospin setup described above develops an instability at large 
enough isospin density. This means that the modes X and Y develop quasi-normal 
modes (QNM) which have a positive imaginary part, i. e. which are enhanced in- 
stead of being damped. This suggests that exactly the mesons corresponding to the 
X, Y fluctuations condense. This instability comes about naturally if we think about 
the fact, that we are trying to push more and more charge density into a confined 



Fig. 6 A comparison between 
the finite temperature part of 
the spectral functions Mxy 
and 2%yx (solid lines) in 
the two flavor directions X 
and Y transversal to the chem- 
ical potential is shown in 
units of N c T 2 T r /4 for large 
quark mass to temperature 
ratio xo = 0.99 and d = 0.25. 
The spectral function M E igi 
along the a = 3-flavor direc- 
tion is shown as a dashed 
line. 
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volume. In particular the 3-7 strings charging the D7-brane are located at the black 
hole horizon as motivated earlier in section 13.11 Our current setup does not allow 
the strings to move into the bulk because we are forcing all the background fields 
they would create there to be zero. Up to now we have required Aq ^ and A^ = 
for all other field components. But we are going to relax that restriction by allow- 
ing a non-trivial A3. We will see below that this is sufficient to stabilize the theory 
in a new phase which we will prove to be superconducting/superfluid. This setup 
naturally produces a p-wave structure since we have seen above that the condensing 
vector mesons have a triplet structure. According to tableQ]this implies a p-wave (at 
least to leading order). 



5 Signatures of super-something 

Finally we put together everything we have learned about D-branes, superconduc- 
tors and holographic methods. We discuss the holographic results and interprete 
them in the condensed matter context. Section 15.11 starts with the thermodynam- 
ics and section IBT21 continues with the details of the fluctuation computation. The 
spectrum and conductivity is examined in section IB31 and all the signatures will be 
pointing to the fact that we have indeed created a holographic p-wave superconduc- 
tor/ 'superfluid. 



5.1 Thermodynamics of the broken phase 

FigureQshows the background field configuration. The different curves correspond 
to the temperatures T = T c and T w Q.9T C . The plots are obtained at zero quark 
mass m = and by using the adapted symmetrized trace prescription. Similar plots 
may also be obtained at finite mass m ^ and by using the DBI action expanded to 
fourth order in F. These plots show the same features: (top left) The gauge field Aq 
increases monotonically towards the boundary. At the boundary, its value is given 
by the dimensionless chemical potential fi. (top right) The gauge field A3 is zero for 
T > T c . For T < T c , its value is non-zero at the horizon and decreases monotonically 
towards the boundary where its value has to be zero, (bottom left) The conjugate 
momentum /5q of the gauge field Aq is constant for T > T c . For T < T c , its value 
increases monotonically towards the boundary. Its boundary value is given by the 
dimensionless density d^. (bottom right) The conjugate momentum p\ of the gauge 
field A3 is zero for T > T c . For T < T c , its value increases monotonically towards 
the boundary. Its boundary value is given by the dimensionless density — d\ . 

All thermodynamic quantities are determined in terms of the relevant thermody- 
namic potential. We use the grand potential W7 in the grandcanonical and the free 
energy F-j in the canonical ensemble. Both stem from the D7-brane action. They are 
related through a Legendre transformation in the background gauge field Aq on the 



28 



Matthias Kaminski 



gravity side. We obtain both potentials due to the gauge/gravity dictionary from the 
Euclideanized gravity on-shell action according to Z = e~ /on - shl:11 , with the partition 
function Z of the boundary field theory. So the for example W 7 = TS^flf^f. 
Comparison of the grand potentials in figure[8]then shows that the new phase with a 
finite value for A3 is thermodynamically preferred below a temperature T = T c and 
does not exist above. The transition seems numerically smooth, i. e. it is a continu- 
ous phase transition. In these respects both computation schemes, the symmetrized 
trace prescription as well as the DBI expansion to fourth order give the same qual- 
itative behavior. Figure [9] confirms this result and at least numerically determines 
the transition to be of second order. The order parameter d\ vanishes with a criti- 
cal exponent of 1/2 to numerical accuracy. It has been explicitly verified that this 
condensation of vector particles occurs before, i. e. at higher temperature than the 
condensation of the scalars in our theory (see remarks in J6]]). If this was QCD the 
scalar pions would condense before the vector mesons did. However, there might be 
other instabilities and also other symmetry breaking configurations may be possible. 

Using some intuition to define the superconducting density d s = (Jq — cq)/Jq, 
we find that it vanishes linearly at the critical temperature as shown in figure [Tol 

All these signatures are those of a superconducting/superfluid phase transi- 
tion. The order parameter d\ has vector structure by construction, which im- 
plies the p-wave. 




Fig. 7 Profiles of the relevant dimensionless gauge fields A on the D7-branes and their dimension- 
less conjugate momenta p versus the dimensionless AdS radial coordinate p near the horizon at 
p = 1. (top left) the lower curve is T = 0.9T C , (bottom left) the constant is at T = T c since there are 
no charges in the bulk, the other curve shows that the charge proportional to increases towards 
the boundary. 
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We can also compute the specific heat which the flavor branes contribute to the 
theory as Cj = —Td 2 Fq/dT 2 . In figure[TT]the higher curve corresponds to the nor- 
mal phase with A3 = 0, while the lower one corresponds to the superconducting 
phase with A\ ^ 0. Note that the total specific heat is always positive although the 
flavor brane contribution is negative. The divergences near T = in both phases 
can be attributed to the missing backreaction in our setups We read off from the 
numerical result that near the critical temperature, the dimensionless specific heat 
%= 32Ct /(XNfN c T 3 ) is constant in the superconducting phase. 

This implies that the dimensionful specific heat C7 is proportional to T 3 . This 
temperature dependence is characteristic for Bose liquids. There is also a finite 
jump in the specific heat at T c . 



T/T c T/T c 

Fig. 8 Grand canonical potential computed (a) from the symmetrized trace prescription and (b) 
from the expanded DBI action, both at vanishing quark mass M q — 0. The qualitative behavior 
agrees for both prescriptions: At T = T c the energy curve splits into a lower and a higher energy 
branch. The branch with lower energy in both cases is the one with a finite new condensate d\ 
below T c . 

Fig. 9 The order parameter 
d\ defined in 452t as obtained E0 
from the adapted symmetrized 
trace prescription versus 50 
temperature: The case of 
vanishing quark mass (top 
solid curve) shows the same 30 
behavior near T c as that where 
fi/Mq = 3 is fixed (lower solid 20 
curve). Both curves go to zero 
with a critical exponent of 1/2 
near T c , as visualized by the 
fit 55(1 -T/T L .) 1 / 2 (dashed 
curve). 




6 The backreaction of the gauge field on the geometry, i. e. on the Einstein equations for metric 
components has been considered in 1 49 1 . 
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5.2 Fluctuations in the broken phase 



Let us now investigate fluctuations of this system. We are going to see that the 
computation of the non-Abelian DBI-action is a very subtle issue. We argue for a 
novel prescription which actually gives reasonable results. The full gauge field A on 
the branes consists of the field A and fluctuations a, 



-A^df- 



-A\t 1 &ci- 



(73) 



where x a are the St/ (2) generators. The linearized equations of motion for the fluc- 
tuations a are obtained by expanding the DBI action in a to second order. We will 
analyze the fluctuations a\ and X = a\ + \a 



Y = a\ — ia\. 



Including these fluctuations, the DBI action reads 



S=-T 7 J d 8 £ Str y' det [G + (2na')F] , 
with the non-Abelian field strength tensor 

pa =F a ,pa 

where the background is collected in 



(74) 



(75) 



Fig. 10 Superconducting 
density d s = — co)/Jq 
versus temperature T: In 
both, the massless and the 
massive case at [l/M q = 3 
(both curves coincide on the 
plot), the superconducting 
density d s vanishes linearly 
at the critical temperature. 
This is visualized by the fit 
6.8(1 -T/T c ) (dashed). 



Fig. 11 The flavor brane 
contribution to the specific 
heat as computed from the 
adapted symmetrized trace 
prescription in the massless 
case. The lower curve starting 
at T = T c is the heat capacity 
in the superconducting phase. 
The other curve is computed 
in the normal phase. 
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]7 a —12, A" -l J—f abc A b A c 

t HV — la [n A v] + jjJ A |i v ' 

and all terms containing fluctuations in the gauge field are summed in 

= 2<Vv] + -j%f abCa >v + -^f abc K a v + 4 A v) • ( 77 > 

Index anti-symmetrization is always defined with a factor of two in the following 

way <5[ M A V] = (d^A v - d v A^)/2. 



5.2.1 Adapted symmetrized trace prescription 

In this section we use the adapted symmetrized trace prescription to determine the 
fluctuations about the background we discussed in section U.3.11 To obtain the lin- 
earized equations of motion for the fluctuations a, we expand the action (|74| | to 
second order in fluctuations, 



S (2) 



-T 7 J d 8 £Str 
(2na') 2 r 



{2na') 



SfSf"" F^yi 



{2na') 2 



(78) 



V'jX 



As in |42J, we collect the metric and gauge field background in the tensor Sf = 
G + (2na')F. Using the Euler-Lagrange equation, we get the linearized equation of 
motion for fluctuations at in the form 



= d A Str 
-Sir 



1 



_^d^[ K X] + ( 2w oO(sf M[r Sf A]v + ^ y & [Kl] )F Vfl } 
f bd x"V^x {sf (a +A)'{ + {2na') (V [ ^ ]v 



(79) 



Note that the linearized version of the fluctuation field strength used in equation d79l i 
is given by 



F" = 2d u 



* rabc { a b 



^a a v] + -j=f {A^ay + a»^ v ) + 0{a L ) . (80) 
In our specific case the background tensor in its covariant form is given by 

%v = G^ v t° + (2na')(2dpAl8 4 ^8 v] oT 3 +2dpA\8 A[tl 8 v] iZ l +2-^4^<5 0[m <S v]3 t 2 

(81) 
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Inversion yields the contravariant form needed to compute the explicit equations of 
motion. The inverse of Sf is defined as S^ v Sf VM ' = S£,T°|H The non-zero compo- 
nents of ^ v may be found in the appendix of |6j. 

Fluctuations in a\: For the fluctuation a\ with zero spatial momentum, we obtain 
the equation of motion 

o = («!)"+ ¥(«!)'- 



4p» , //h2 



16- 



P 4 /- 



HI 



2c 2 



(82) 



3 = JA\/{2V2%) and H = V^G 22 ^ 44 



with , 



Fluctuations in X = a\ 



ia^: For the fluctuations X and Y with zero 



spatial momentum, we obtain the coupled equations of motion 



d2 d 



6f 33 
2 f ^ 44 



«{03} 



«4 



2p^p[v^G 22 ^ 04 >(v 1 >+^ 3 )1 



R 2 



.#3 ^# Q 3 



(83) 



where the component of the inverse background tensor may be found in the ap- 
pendix of |6 1 (just like the corresponding formula for Y which is only different from 
( 1831 ) by a few signs), index symmetrization is defined SfW — (cf'J + <^j')/2 and 
J*l = rAl/(2V2n). 



5.2.2 Expansion of the DBI action 

In this section we determine the equation of motion for the fluctuation a\ in the 
background determined by the DBI action expanded to fourth order in F (see section 
13.3.2b . To obtain the quadratic action in the field a|> we fi rst have to expand the DBI 
action d74l i to fourth order in the full gauge field strength F, and expand the result to 
second order in a\. Due to the symmetries of our setup, the equation of motion for 
the fluctuation a\ at zero spatial momentum decouples from the other equations of 
motion, such that we can write down an effective Lagrangian for the fluctuation a 2 . 
This effective Lagrangian is given in the appendix of (6). The equation of motion 



7 We calculate the inverse of 'S by ignoring the commutation relation of the t's because of the sym- 
metrized trace. It is important that z"x b must not be simplified to e abc X c since the symmetrization 
is not the same for these two expressions. 
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for a\ with zero spatial momentum determined by the Euler-Lagrange equation is 
given by 



o=(4)"+ d - E f(4)'-§ 



4 ^ w 




8 d p [^GG m G 22 G^G^Al{d p Al)(^ 

"3 3? 



'2 • 



(84) 
(85) 



where 3tf = ^/^GG 22 J? 44 . We introduce the factors Jf'j which may be found in 
the appendix of Q to emphasize the similarity to the equation of motion obtained 
by the adapted symmetrized trace prescription i 



5.3 Conductivity & spectrum 



We calculate the frequency-dependent conductivity a(a>) using the Kubo formula, 



a((0) = —G R ((O,q = 0), 

CO 



(86) 



where G R is the retarded Green function of the current dual to the fluctuation 
a\, which we calculate using the method obtained in P4l . The current j\ is the 
analog to the electric current since it is charged under the 1/(1)3 symmetry. In real 
space it is transverse to the condensate. Since this fluctuation is the only one which 
transforms as a vector under the 50(2) rotational symmetry, it decouples from the 
other fluctuations of the system. 

Exercise: Prove equation A86[ for the current j\ assuming that the gauge/gravity 
correspondence is correct. Recall that in regular electrodynamics a = J/E with 
the electric current J and the electric field E. Also recall that the two-point Green 
function for a current J dual to the gauge field A can holographically be written in 
the form G R °c d p A/A. 




Fig. 12 Conductivity (a) from the symmetrized trace prescription at T/T c = 10, 1, 0.5, 0.28 (from 
left to right), (b) from the DBI action expanded to fourth order at T/T c — °°, 1, 0.6, 0.5, 0.39. 
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(Pseudo) gap An analysis of the imaginary part of the conductivity using Kramers- 
Kronig realtions shows a delta peak at vanishing frequency w = in the real part. 
The frequency-dependent conductivity is shown for our two distinct computation 
schemes in figure [12] Independent from the scheme we see an energy gap develop 
and grow while the temperature is decreased. The temperature of the black hole hori- 
zon induced on the D7-branes is proportional to the inverse particle mass parameter 
rrT 1 ex; T. Therefore from the trivial flat brane embedding at m = we get an infinite 
temperature in subfigure[T2](b). Both schemes show the development of peaks in the 
conductivities. The peaks coming from fluctuations around the symmetrized trace 
prescription background are a lot more pronounced. Taking into account only sec- 
ond order terms in the expanded DBI action as in [40 1 would hide the peak structures 
completely. Therefore we conclude that these peaks are higher order effects. Since 
the conductivity is closely related to the spectral functions, we interprete the peaks 
as quasiparticles just as described in section |4. 1 .21 and l4~2l In particular these are 
again vector mesons. This identification is confirmed by the spectral function in fig- 
ureQj] The peaks in that figure are identical with those in the conductivity and they 
approach the supersymmetric line spectrum for vector mesons, as in section l4~2l 

Dynamical mass generation Even if we choose the two D7 branes to coincide with 
the stack of D3-branes, i. e. if we choose the quark mass to vanish, we observe the 
quasi-particle peaks mentioned before. This is due to a Higgs-like mechanism which 
dynamically generates masses for the bulk field fluctuations, which in turn give 
massive quasi-particles in the boundary theory. In the bulk our fields Aq and A3 break 
the SU(2) symmetry spontaneously since the bulk action is still SU (2) -invariant. 
Thus there are three Nambu-Goldstone bosons which are immediately eaten by the 
bulk gauge fields, giving them mass. This can be seen explicitly in the action where 
the following mass terms for the gauge field fluctuations appear:(Ao) 2 (a 1,2 ) 2 and 
(A')V) 2 . 



Fig. 13 Finite temperature 
part of the spectral function 
- M Q with ,'M = Altwhn 
units of NfN c T 2 /& versus 
the dimensionless frequency 
w = (o/(27tT) at finite quark 
mass m = 2.842 and chemical 
potential pL = 3.483. The 
grey lines correspond to 
the supersymmetric mass 
spectrum calculated in 1241 . 
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The Meissner effect is a distinct signature of conventional and unconventional su- 
perconductors. It is the phenomenon of expulsion of external magnetic fields. An 
induced current in the superconductor generates a magnetic field counter-acting the 
external magnetic field H. In AdS/CFT we are not able to observe the generation 
of counter-fields since the symmetries on the boundary are always global. Never- 
theless, we can study their cause, i.e. the current induced in the superconductor. As 
usual ||50l |5T] [52 9| the philosophy here is to weakly gauge the boundary theory 
afterwards. In order to investigate how an external magnetic field influences our p- 
wave superconductor, we have two choices. Either we introduce the field along the 
spatial z-direction 7/|t 3 or equivalently one along the y-direction i.e. //|t 3 . Both 
are "aligned" with the spontaneously broken t/(l)3-flavor direction. 

As an example here we choose a non-vanishing // 3 T 3 . This requires inclusion 
of some more non-vanishing field strength components in addition to those given in 
equation d26l). In particular we choose A^(p ,x) , A\{p ,x) andA 3 (x) =x l yielding 
the additional components^] 

Hl=Fl 2 = -Fl = diA\, (87) 

= = ^4*3. (89) 
*h> = ~Fl = diAl (90) 

Recall that the radial AdS-direction is designated by the indices p or 4 synony- 
mously. Amending the DBI-action ( f3Tb with the additional components d87l l. we 
compute the determinant in analogy to equation ( T3Tb . We then choose to expand the 
new action to second order in F, i.e. we only consider terms being at most quadratic 



Fig. 14 The line of critical 
magnetic field versus critical 
temperature. Below this line 
the external magnetic field 
coexists with the supercon- 
ducting condensate. Above 
the line the superconducting 
condensate vanishes. We set 
the spatial position arbitrarily 
to x = 0. 1 since the critical 
line does not depend on x. 



11 



T 



Close to the phase transition, it is consistent to drop the dependence of the field Hj on p. Away 
from the phase transition the p dependence must be included. From the boundary asymptotics it 
will be possible to extract the magnetic field and the magnetization of the superconductor. 
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in the fields. This procedure gives the truncated DBI action 



SDBi = -T m Nffd*ZV=G 



1 



+G 00 G 44 (F p 3 ) 2 + (G»r(Ft 3 y + (G»Y(F&) 



33\2/ r l \2 



G 00 G 33 (F 1 3 ) 2 + (G-)^(^ 2 ) 



33N2/ r 3 \2 



(91) 



Respecting the symmetries and variable dependencies in our specific system, this 
can be written as 



Sdbi = ~Tw N f / d 8 £ y/=G 



1 



1 + -(G 00 G 33 (^A 3 ) 2 + G^G 44 (^) 



(92) 



+G 00 G 44 («3 p A^ + (G 33 ) Wa)* + (G s r(H() 



J3\2 



33\2/3 Tl\2 



-33n2/£>3\2 



r 00 r 33 PH T (pXl\2 

+G G (2^X (AoA3) 



with the convenient redefinitions 



4„2 n 4„,2 

33 s2 PH 7 



(G 3 



(2na') 2 X 



2nd A _, , , 

A= A, x = p H x, H 3 =27tdH 3 , p = pnP- 

Ph 



(93) 



Rescaling the x-coordinate once more 



„ = mpH 2 j - 

the equations of motion derived from the action d92l take a simple form 

a _ ,3221 , 3 p [^GG^] d ri r-Hf (£3 r,2Tl _ =2 711 _ J^G™ / T3^2T 



(94) 



(95) 



Here all metric components are to be evaluated at R = 1 and p — > p . 

We aim at decoupling and solving the system of partial differential equations 
by the product ansatz 



A](p,x) = v(p) M (i). 



(96) 



For this ansatz to work, we need to make two assumptions: First we assume that 
Ag is constant in x. Second we assume that A3 is small, which clearly is the case 
near the transition T — > T c . Our second assumption prevents Aq from receiving a 
dependence on.f through its coupling to A\(p,x). These assumptions allow to write 
the second equation in ((95) as 
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J p"\V> ' v^GG^G 44 "P"^' 27? G 4 

-^^v(p)^^. (97) 



All terms but the last one are independent of x, so the product ansatz 
tent only if 

f»W =c , (98) 
u{x) 

where C is a constant. The differential equation d98l has a particular solution if 
C = — (2n + 1) ,« € N. The solutions for w(Jc) are Hermite functions 

«„(*) = ! H n (x),H n (x) = (-l)" e %iV^, (99) 

yn!2"\/7r " x 

which have Gaussian decay at large 3> 1. Choosing the lowest solution with n = 
and Hq = 1, which has no nodes, is most likely to give the configuration with lowest 
energy content. So the system we need to solve is finally given by 



B p [V=GG°°G«] , _ 3 f G 33 



= ^Ag+ ^ G oo G 44 J ^o^2^o(«o(^)v(p)r, (100) 
«9„ [\/^GG 33 G 44 l , y 2 G 00 ,~,, 9 y// 3 G 33 

Asymptotically near the horizon the fields take the form 

Al = + fl2 (p-l) 2 + ^((p-l) 3 ), (102) 
v=b +^(p-l) 2 + ff((p-l)- 3 ), (103) 
while at the boundary we obtain 

Al=fi + f| + ^(p- 4 ), (104) 
P 

v = + ^ + 6\p- A ). (105) 
P l 

We succeed in finding numerical solutions v(p) and Ajj(p) to the set of equations 
( llOOl l obeying the asymptotics given by equations ( 1102b and ( 11041 i. These numer- 
ical solutions are used to approach the phase transition from the superconducting 
phase by increasing the magnetic field. We map out the line of critical temperature- 
magnetic field pairs in figure [14] In this way we obtain a phase diagram displaying 
the Meissner effect. The critical line in figure [14] separates the phase with and with- 
out superconducting condensate d\ . 
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We emphasize that this is a background calculation involving no fluctuations. 
Complementary to the procedure described above we also confirmed the phase dia- 
gram using the instability of the normal phase against fluctuations. Starting at large 
magnetic field and vanishing condensate d\, we determine for a given magnetic 
field //| the temperature T C (H%) at which the fluctuation a\ becomes unstable. That 
instability signals the condensation process into the superconducting phase. 

The presence of the coexistence phase below the critical line, where the system is 
still superconducting despite the presence of an external magnetic field, is the signal 
of the Meissner effect in the case of a global symmetry considered here. If we now 
weakly gauged the flavor symmetry at the boundary, the superconducting current 

would generate a magnetic field opposite to the external field. Thus the phase 
observed is a necessary condition in the case of a global symmetry for finding the 
standard Meissner effect when gauging the symmetry. 



6 Interpretation & Conclusion 

Our findings merge to a string theory picture of the pairing mechanism and the 
subsequent condensation process in section 16.11 Finally we summarize what we 
have learned about holographic p-wave super-somethings and propose some future 
territories to be conquered. 



6.1 String Theory Picture 

We now develop a string theory interpretation, i. e. a geometrical picture, of the 
formation of the superconducting/superfluid phase, for which the field theory is dis- 
cussed in section I2TT1 We show that the system is stabilized by dynamically gen- 
erating a non-zero vev of the current component j\ dual to the gauge field A3 on 
the brane. Moreover, we find a geometrical picture of the pairing mechanism which 
forms the condensate (/J), the Cooper pairs. Let us first describe the unstable con- 
figuration in absence of the field A3. As known from ll37l[3Tll20l . the non-zero field 
Ag is induced by fundamental strings which are stretched from the D7-brane to the 
horizon of the black hole. In the subsequent we call these strings 'horizon strings'. 
Since the tension of these strings would increase as they move to the boundary, 
they are localized at the horizon, i. e. the horizon is effectively charged under the 
isospin charge given by (fl~|i. By increasing the horizon string density, the isospin 
charge on the D7-brane at the horizon and therefore the energy of the system grows. 
In ll20l . the critical density was found beyond which this setup becomes unsta- 
ble. In this case, the strings would prefer to move towards the boundary due to 
the repulsive force on their charged endpoints generated by the flavorelectric field 
4=% = -dpA%. 
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The setup is now stabilized by the new non-zero field A3. We can think of this 
field as being induced by D7-D7 strings moving in the x 3 direction. This movement 
of the strings may be interpreted as a current in x 3 direction which induces the 
magnetic field B\ p = F3L = —d p A\. Moreover, the non-Abelian interaction between 
the D7-D7 strings and the horizon strings induces a flavorelectric field E\ = F 3 2 = 



From the profile of the gauge fields and their conjugate momenta (see figure |7|i 
we obtain the following: For A3 = 0, i. e. in the normal phase (T > T c ), the isospin 
density is exclusively generated at the horizon by the horizon strings. This can 
also be understood by the profile of the conjugate momenta (see figure [7] (bot- 
tom left)). We interpret p^{p*) as the isospin charge located between the horizon 
at p = 1 and a fictitious boundary at p = p*. In the normal phase, the momentum 
Pq is constant along the radial direction p (see figure |7Jbottom left)), and therefore 
the isospin density is exclusively generated at the horizon. In the superconducting 
phase where A3 7^ 0, i. e. T < T c , the momentum p$ is not constant any more. Its 
value increases monotonically towards the boundary and asymptotes to (see fig- 
ure [TJbottom left)). Thus the isospin charge is also generated in the bulk and not 
only at the horizon. This decreases the isospin charge at the horizon and stabilizes 
the system. 

Now we describe the string dynamics which distributes the isospin charge into 
the bulk. Since the field A3 induced by the D7-D7 strings is non-zero in the super- 
conducting phase, these strings must be responsible for stabilizing this phase. In the 
normal phase, there are only horizon strings. In the superconducting phase, some 
of these strings recombine to form D7-D7 strings which correspond to the non-zero 



Fig. 15 Sketch of our string setup: The figure shows the two coincident D7 branes stretched from 
the black hole horizon to the boundary as a green and a blue plane, respectively. Strings spanned 
from the horizon of the AdS black hole to the D7-branes induce a charge at the horizon 1 20 3T1I37I . 
However, above a critical charge density, the strings charging the horizon recombine to D7-D7 
strings. These D7-D7 strings are shown in the figure. Whereas the fundamental strings stretched 
between the horizon and the D7-brane are localized near the horizon, the D7-D7 strings propagate 
into the bulk balancing the flavorelectric and gravitational, i. e. tension forces (see text). Thus 
these D7-D7 strings distribute the isospin charges along the AdS radial coordinate, leading to 
a stable configuration of reduced energy. This configuration of D7-D7 strings corresponds to a 
superconducting condensate. 
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gauge field A3 and carry isospin charge^. There are two forces acting on the D7-D7 
strings, the flavorelectric force induced by the field Ep and the gravitational force 
between the strings and the black hole. The flavorelectric force points to the bound- 
ary while the gravitational force points to the horizon. The gravitational force is 
determined by the change in effective string tension, which contains the p depen- 
dent warp factor. The position of the D7-D7 strings is determined by the equilibrium 
of these two forces. Therefore the D7-D7 strings propagate from the horizon into 
the bulk and distribute the isospin charge. 

Since the D7-D7 strings induce the field A3, they also generate the density d\ 
dual to the condensate (j\ [), the Cooper pairs. This density d\ is proportional to 
the D7-D7 strings located in the bulk, in the same way as the density counts 
the strings which carry isospin charge 1311 . This suggest that we can also interpret 
p\(p*) as the number of D7-D7 strings which are located between the horizon at 
p = 1 and the fictitious boundary at p = p* . The momentum p\ is always zero at the 
horizon and increases monotonically in the bulk (see figure [7] (bottom right)). Thus 
there are no D7-D7 strings at the horizon, nevertheless most of them are located 
near the horizon. 

The double importance of the D7-D7 strings is given by the fact that they are 
both responsible for stabilizing the superconducting phase by lowering the 
isospin charge density at the horizon, as well as being the dual of the Cooper 
pairs since they break the f/( 1)3 symmetry. In QCD-language they correspond 
to quarks pairing up to form charged vector mesons which condense subse- 
quently. 



6.2 Summary 

In conclusion we have derived from the top (string theory) down to the gravity the- 
ory a holographic p-wave superconductor. Thus we were able to directly identify 
the degrees of freedom in the boundary field theory, allowing us to translate ge- 
ometric features directly into field theory features. In particular we have found a 
string theory picture of the pairing mechanism. The Cooper pairs are modeled by 
strings spanned between the two flavor D7-branes corresponding to quasi-particles 
in the vector bi-fundamental representation, i. e. vector mesons. The dual thermal 
field theory is 3+1 -dimensional JV = 2 supersymmetric Yang-Mills theory with 
SU (N c ) color and SU (2) flavor symmetry coupled to an JV = 4 gauge multiplet. It 
shows a conductivity gap at low temperatures. A pseudo-gap forms even above T c . 
The onset of the Meissner-Ochsenfeld effect is visible and in the conductivity spec- 



9 Note that the D7-D7 strings are of the same order as the horizon strings, namely Nf/N c , since 
they originate from the DBI action |37| . 
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tram we find massive quasi-particles even at vanishing quark mass. Their masses 
are generated through a Higgs-like mechanism in the bulk. 

Note that our results can also be interpreted using this very setup as a model for 
the quark gluon plasma as introduced in section 13.11 In that case we have found a 
flavor superfluid phase. This should not be confused with the color-superconducting 
phase theoretically found in QCD at high baryon density. 



6.3 Outlook 

Some directions for promising future investigation are the study of critical expo- 
nents near the phase transition. Transport coefficients such as the speeds of second, 
fourth and other sounds can be determined. For some of these it will be neccessary 
or convenient to use the backreacted solutions for this D3/D7 system given in ll49ll . 
Similar to earlier semi-classical drag computations it might be instructive to com- 
pute the drag on the various strings near the superconducting phase transition. This 
computation had been suggested already in [5 1 and carried out for a scalar conden- 
sate superfluid in |53l . It would be interesting to see what happens to the Fermi 
surface formed by the background fermions when the superconducting/superfluid 
phase is entered. Adventurous spirits may also take this setup as a serious model 
for p-wave superconductors such as the ruthenate compounds mentioned in the in- 
troduction. The geometric insights we gain from the dual gravity description can in 
principle be directly translated to precise field theory statements. Again this pleas- 
ant fact is due to knowing the field theory degrees of freedom exactly by using the 
gauge/gravity correspondence. A different setup in which vector mesons condense 
at finite isospin potential is the Sakai-Sugimoto model as shown in l54l . It could be 
illuminating to study all aforementioned effects in this model as well ([55] studies a 
possibly superconducting phase in the Sakai-Sugimoto model but considers a pion 
condensate at finite baryon density). 
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